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Standard  Form  298  (Rev.  2-89) 
Prescribed  by  ANSI  Std  239-18 
298-102 


METHODS  FOR  CONDITIONING  ELECTRON  BEAMS 
IN  FREE  ELECTRON  LASERS 


Many  coherent  radiation  generation  mechanisms  are  based  on  the 
longitudinal  bunching  of  electron  beams.  These  sources  include  traveling 
wave  tubes  and  free  electron  lasers  (FELs)  [1-7].  The  degree  to  which  an 
electron  beam  can  be  bunched  is  a  strong  function  of  the  beam  quality.  The 
two  independent  contributions  to  the  electron  beam  quality  are  the 
intrinsic  energy  spread  and  emittance,  both  of  which  lead  to  a  spread  in 
the  axial  electron  velocity  and  limit  the  operating  wavelength,  gain  and 
efficiency  of  the  device  [7-10].  A  method  for  conditioning,  i.e.,  reducing 
the  axial  beam  velocity  spread,  was  recently  proposed  in  which  the  beam  was 
propagated  through  a  periodic  array  of  focusing,  drift,  defocusing,  drift 
channels  and  microwave  cavities  excited  in  the  non-axially  symmetric  TM21q 
mode  [11]. 

In  this  letter  we  propose  an  alternative  conditioning  method  which 
redistributes  the  electrons'  energy  according  to  their  betatron  amplitude 
by  using  the  electric  field  of  an  axially  symmetric,  slow,  TM  waveguide 
mode.  In  addition,  a  second  method  is  briefly  described  in  which  the 
energy  redistribution  is  correlated  with  the  electrons'  synchrotron 
amplitude. 

We  first  obtain  the  restriction  on  the  operating  FEL  wavelength  due  to 

emittance  and  energy  spread.  In  the  FEL  mechanism,  the  resonance  condition 

isw-v  (k  +  k  )  =  0,  where  w  =  ck  is  the  frequency,  v  is  the  axial 
z  w  z 

electron  velocity,  k  =  2it/X  and  X  is  the  wiggler  wavelength.  This 

w  w  w 

condition  implies  that  the  beam's  axial  velocity  spread,  Sv  ,  should 
Manuscript  approved  January  15,  1993. 
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satisfy  S3  «  X/2L  where  S3,  =  Sv  /c,  X  is  the  radiation  wavelength  and  L 
J  z  z  z 

is  the  interaction  length  (e-folding  length)  of  the  radiation  field  in  the 
low  gain  (high  gain)  regime.  The  axial  velocity  spread  can  be  written  as 

S3z  =  [(1  +  a^/2)Sy/y  -  e2/2r2)/y2  [see  Eq.  (2)],  where  y  -  1  +  E/mQc2  is 
the  relativistic  factor,  E  is  the  beam  energy,  Sy/y  is  the  fractional 
intrinsic  beam  energy  spread,  en  is  the  normalized  emittance  [12],  ay  is 
the  wiggler  strength  parameter  and  rb  is  the  radius  of  the  matched  electron 
beam.  Using  this  expression  for  S3Z  we  find  the  restriction  on  the 
operating  wavelength, 


X  »  (ji)1/2(L/Zr)" 


1/2 

e2/2r?  -  ( 

1  +  a2/2)Sy/Y 

Y 

n  b  \ 

k  w  ) 

1/2 


(1) 


where  Z  =  Jtr2/X  is  the  vacuum  Rayleigh  length  associated  with  a  Gaussian 
R  s 

2  2. 

beam  with  minimum  spot  size  rg  and  f  =  rb/rg  is  the  filling  factor. 

The  inequality  in  (1)  can  be  simplified  by  noting  that,  in  both  the 
low  and  high  gain  regimes  of  the  FEL,  L  =  ZR  for  f  =  1  [13-16].  That  is, 
in  the  low  gain  regime,  L  =  ZR  and  f  =  1  are  required  to  minimize 
diffraction  effects  and  maximize  gain.  In  the  high  gain  optically-guided 
regime,  it  can  be  shown  that  the  e-folding  length  is  approximately  equal  to 
the  vacuum  Rayleigh  length,  i.e.,  L  =  ZR,  for  f  -  1.  In  many  cases, 
electron  beam  quality  is  limited  by  the  emittance  contribution  and  not  the 
energy  spread  term,  i.e.,  6y/y  <<  (l/2)(en/rb)  .  In  this  case,  the 
wavelength  limit  in  FELs  can  be  simply  written  as  X  »  e^y.  It  is  clear 
that  electron  beam  quality,  in  particular,  &3z,  limits  the  operating 
wavelength  of  FELs. 

To  analyze  our  conditioning  methods  we  first  consider  the  electron 
trajectories  in  a  planar  wiggler  with  parabolic  pole  faces.  These  orbits 
consist  of  rapidly  varying  (wiggler  period  scale  length)  and  slowly  varying 
(betatron  period  scale  length)  terms  [17].  For  the  wiggler  field 
components  used  in  [17],  the  approximate  electron  orbits  are  x  =  xf  +  xg 


N 
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and  y  =  y.  +  y  .  The  fast  components  of  the  orbit  are  x^(z)  -  (av/rhw>(l 

i  ®  2 

+  k2(x2  +  y2)/4)cos  (k  z),  yf(z)  «  0,  where  ay  =  |e|  Bw/(kwmoc  ),  ky  = 

W  S  S  w  i- 

2n/X  ,  X  is  the  viggler  period,  while  the  slow  components  are  xg(z)  = 
w  w 

xocos  (Kpz  +  0OX),  ys(z)  =  yocos  (KpZ  +  0oy),  where  (xQ,  yQ)  are  the 
maximum  amplitude  of  the  betatron  oscillations  in  the  x  and  y  directions, 
respectively,  K$  -  <k2  -  k2)1'2,  k0  .  avku/2T  «  ky  is  the  betatron 

wavenumber,  k2  -  (m2/c2)/(2t3>,  a2  =  4kM2VV  no  is  thc  b€aI"  de"Sity 

and  G  ,  G  are  constant  phases.  Ve  have  assumed  that  aw/Y  «  ^wrb  <<  1 
ox  oy 

and  rk  rL  «  1  which  implies  that  v/Y  «  1  where  v  =  I[A]/17000Pz  = 

P  b 

W2r2/4C2  Budker's  parameter  and  I  is  the  current  in  amperes, 
p  b 

In  the  highly  relativistic  limit,  the  axial  particle  velocity 

2  2  2 

normalized  to  the  speed  of  light  is  given  by  0Z  =  1  -  1/2y  -  Ox  +  Py)/2‘ 

As  pointed  out  in  Ref.  17,  the  square  of  the  perpendicular  velocity, 
averaged  over  the  wiggler  period,  is  independent  of  z.  Substituting  the 
fast  and  slow  orbits  into  the  expression  for  Pz  and  setting  Y  =  rQ  +  &Y, 
where  5y  is  the  electron's  energy  deviation  term,  we  find  that  Pz  =  PQZ  + 
SPz,  with  eoz  =  i  -  <1  +  av/2)/<2^o>»  and 

SPZ  -  (l  ♦  a2/2)sy/T3  -  k2  r2/2.  (2) 

Here  y  is  the  gamma  associated  with  the  reference  electron,  traveling 
along  the  z-axis  without  a  betatron  oscillation,  r2  =  xq  +  yQ,  and  terms 
varying  on  the  wiggler  wavelength  scale  have  been  neglected.  The  space 

O 

charge  contribution  to  Pz  is  of  order  v/yq  and  is  neglected.  The 
normalized  beam  emittance,  for  a  matched  beam  in  the  focusing  fields  of  the 
wiggler,  is  en  =  Y^r2.  Note  that  the  emittance  contribution  to  the 
velocity  spread  in  (2),  i.e.,  k2r2/2  is  independent  of  propagation 
distance.  It  will  be  assumed  that  the  axial  velocity  spread  due  to 
emittance  initially  dominates  the  velocity  spread  caused  by  the  intrinsic 
energy  spread. 
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The  proposed  conditioning  field  is  an  axially  symmetric,  slow,  TM 
waveguide  mode  (Fig.  1)  with  axial  electric  field 

E  =  -  E  I  (k,  r)  cos  <f/,  (3) 

z  o  o  1 


together  with  the  associated  transverse  electric  and  magnetic  fields,  where 

E  is  the  maximum  electric  field  amplitude  on  axis,  k.  is  the  transverse 
o  ± 

2  2  1/2 

wavenumber,  k  is  the  axial  wavenumber,  <i)  =  c(k  -  k^)  is  the  frequency, 
i|/  =  kz  -  cot  is  the  phase  and  In  is  the  modified  Bessel  function  of  order  n. 
The  axial  phase  velocity  of  the  traveling  wave  is  matched  to  the  axial  beam 
velocity,  (3  ^  =  w/ck  =  3o,  where  Bo=(1-1/y0)  is  the  normalized  axial 
velocity  of  the  reference  electron.  To  maintain  synchronism  between  the 

conditioning  field  and  the  electrons,  the  axial  and  transverse  wavenumbers 

9-1/2  2  -1/2 

must  satisfy  k  =  (<o/c)ro(r0  -  1)  =  (<a/c)/Bo  an^  k^  =  («/c)(yo  -  1) 

=  (co/c)/(yo0o)»  respectively. 

In  our  conditioning  methods,  the  beam  electrons  are  given  an  energy 
increment  which  cancels  out  the  emittance  contribution  to  the  axial 
velocity  spread.  To  reduce  the  velocity  spread  to  zero,  the  conditioning 
field  must  give  all  the  individual  electrons  a  different  fractional  energy 


increment,  8yc/yq»  given  by 


2.2  2 
Y  kar 
o  P  o 


2  2 
e  r 
n  o 


(i  ♦  4,2) "  2(1  ♦ 


a2/2]ru 
w  ;  b 


(4) 


In  the  first  method,  the  energy  increment  is  proportional  to  the 
square  of  the  betatron  amplitude  and  the  electron  pulse  length  remains 
approximately  constant.  Our  results  show  that  the  degree  of  beam 
conditioning  in  the  first  method  can  be  significantly  improved  by  removing 
the  accelerating  component  of  the  TM  field.  A  straightforward  way  to 
achieve  this  is  to  introduce  an  inductively  generated  axial  electric  field 
in  addition  to  the  waveguide  mode.  This  uniform  inductive  field  may  be 
generated  by  an  azimuthal  magnetic  field  that  is  confined  to  the  interior 
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of  a  cylindrical  shell  of  high-permeability  and  high-resistivity  material, 
e.g.,  a  ferrite. 

The  rate  of  change  of  electron  energy  and  phase  with  respect  to  z  is, 
respectively, 


Y'  =  aQ  ®  ^Io(k^r)cos  -  1  +  |-  r'  I^k^r)  sin  ij/j ,  (5) 

r  =  (a/c)  SB  /p2,  (6) 

z  o 


where  the  prime  denotes  a  derivative  with  respect  to  z  and  aQ  = 

le lEo/(mocw) * 

Conditioning  requires  that  the  spreading  of  the  electron  phase 

(expansion  of  the  prebunched  electron  beam)  be  limited.  An  upper  bound  on 

the  phase  change  is  obtained  by  integrating  Eq.  (6)  with  only  the  emittance 

term  in  &BZ>  l'l'(z)  -  V>(0)  |  <  kk^z/2.  If  the  phase  remains  small,  i.e., 

|<J/|  «  kr/(/2yo)  and  |«|»|  «  kr/(2r^r'),  Eq.  (5)  simplifies  to  y'  = 

ao(<o/c)(kr/2ro)2,  for  k±r  «  1  and  *  «  0.  The  transverse  dynamics  of  the 

electrons  is  governed  by  only  the  wiggler  field  and  the  conditioning  field 

can  be  neglected  if  k2  »  (aQk2/(2Yg))sin  and  k^  »  aQk/(2Y0)- 

Substituting  the  slowly  varying  orbits  into  the  expression  for  y'  and 

2  2 

integrating,  we  obtain  y  =  yQ  +  $r»  where  Sy  =  ((krQ)  /8YQ)ao(»/c)z. 

Making  use  of  Eq.  (2),  the  effective  fractional  axial  energy  spread  is 


f(l  +  a2/2)k3aoz 


4y: 


(7) 


provided  k^z  =  nn,  n  =  1,2...  From  Eq.  (7)  it  follows  that  complete 
conditioning  of  the  beam  is  achieved  at  k^z  =  nn,  provided  the  normalized 
strength  of  the  waveguide  field  is 


a 

o 


4Y>6/k)3  _  y]  a3 
nn(l+a2/2)  2nn  (1  +  a2/2) 


(8) 
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A  full  scale  particle  simulation  of  an  electron  beam  in  a  wiggler  and 
conditioning  field  has  been  carried  out.  The  fully  relativistic  Lorentz 

force  equations  were  integrated  for  10  particles,  using  the  standard 
leapfrog  algorithm  [18]. 

The  first  beam  conditioning  method  is  illustrated  with  two  examples,  a 

10  MeV  and  a  1  MeV  electron  beam,  see  Table  I.  For  the  10  MeV  example,  the 

axial  velocity  spread  of  the  beam  in  the  conditioning  fields  will  reach  a 

minimum  at  z  =  n/k-  =  375  cm,  beyond  which  it  increases  to  its  original 

P 

value.  To  maintain  the  minimum  spread  the  conditioning  field  is 

adiabatically  removed  at  z  <  n/k^.  Figure  2  shows  the  evolution  of  the 

fractional  axial  energy  spread  for  100  randomly  selected  electrons  as  a 

function  of  distance  along  the  waveguide.  The  convergence  of  the 

trajectories  in  Fig.  2  with  propagation  distance  indicates  that  the  spread 

in  axial  velocity  of  the  electrons  is  significantly  reduced  by  the 

conditioning  field.  Figure  3  shows  the  root  mean  square  (rms)  beam  axial 

energy  spread,  )  »  3-S  a  function  of  distance.  In  this  illustration 

o  z  rms 

the  spread  is  reduced  by  a  factor  of  ~  40.  For  the  1  MeV  example,  the  rms 

spread  in  the  axial  velocity  is  observed  to  be  reduced  by  a  factor  of 

approximately  30.  In  both  examples,  the  required  value  of  the  conditioning 
field  is  in  excellent  agreement  with  the  analytical  prediction  in  Eq.  (8). 
For  a  waveguide  diameter  of  1  cm  the  power  in  the  conditioning  field  is  ~ 

10  MW  and  ~  4  kW  for  examples  1  and  2,  respectively. 

In  the  second  method,  the  energy  redistribution  is  correlated  with  the 
electrons'  synchrotron  amplitude.  This  method  can  be  viewed  as  a  rotation 
of  the  electron  distribution  in  phase  space  ( ^').  The  beam  axial 

velocity  spread  and  length  are  proportional  to  the  spread  in  5i|>'  and  Sij/ 

respectively.  Under  certain  conditions  the  beam  distribution  undergoes  a 
rotation  in  phase  space  such  that  the  axial  velocity  spread  decreases  at 
the  expense  of  its  length.  Here,  the  electrons  interact  with  the  wiggler 
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field  and  the  TM  waveguide  field  in  Eq.  (3)  and  the  electron  pulse  is 

phased  such  that  *  =  n/2  +  Sip,  where  |fi*|  «  1  and  k±r  «  1.  The 

electrons'  synchrotron  oscillations  are  described  by  the  pendulum  equation, 

S*»  +  =  0,  where  Kg  «  (ao/yy2z)1,2<*/c  is  the  synchrotron  wavenumber. 

The  beam's  axial  velocity  spread  reaches  a  minimum  after  propagating  a 

distance  L  =  n/2K  (1/4  of  a  synchrotron  period).  By  this  method,  we  have 
s  s 

obtained  reductions  of  >  10  in  the  rms  axial  velocity  spread  in  unoptimized 
simulation  studies. 

In  conclusion,  in  this  paper  two  methods  are  proposed  for  dramatically 
reducing  the  electron  axial  velocity  spread  in  FELs.  A  reduction  in  axial 
velocity  spread  can  improve  the  gain  and  efficiency  of  FELs  as  well  as 
traveling  wave  tubes.  The  beam  conditioning  field  is  that  of  an  axially 
symmetric,  slow,  TM  waveguide  mode.  In  the  first  method,  a  reduction  in 
the  velocity  spread  by  a  factor  of  40  was  obtained.  In  the  second  method, 
reductions  of  >  10  in  the  velocity  spread  have  been  obtained. 
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Table  I 


.J 


Electron  Beam 
Energy,  E 
Emittance,  en 
RMS  Radius,  r^ 

Initial  Axial  Energy  Spread 
Final  (Min.)  Axial  Energy  Spread 


Example  #1 
10  MeV 

3.5x10  cm-rad 
0.14  cm 
2.3xlO~4 
5.8xl0-6 


Example  #2 
1  MeV 

3.5x10 cm-rad 
0.14  cm 
2.4xl0~4 
7.9xl0“6 


Wiggler 

Strength  Parameter,  ay 

Period,  X 
’  w 

Betatron  Period,  X_ 

P 

Conditioning  Field 
Wavelength,  X 
Strength  Parameter,  aQ 
Electric  Field,  E 

o 

Interaction  Length,  ~  X^/2 


0.175 
3.14  cm 
754  cm 


2  cm 

0.1 

160  kV/cm 
375  cm 


0.175 
3.14  cm 
104  cm 


2  cm 
2xl0~3 
3.2  kV/cm 
53  cm 


Table  I.  Simulation  parameters  for  conditioning  a  10  MeV  (Example  #1)  and 
1  MeV  (Example  #2)  electron  beam. 
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TM  Wave  (Axially  Symmetric) 


Slow  Wave  Structure 


Fig.  1  Schematic  of  the  beam  conditioning  configuration  showing  the 
axially  symmetric,  slow,  TM  wave  conditioning  properly  phased 
electron  pulses.  The  phase  velocity  of  the  TM  wave  is  3ph  and  the 
axial  velocity  of  the  reference  electron  is  (both  normalized  to 
the  speed  of  light).  An  inductively  generated  electric  field  may 
be  employed  to  cancel  the  accelerating  component  of  the  TM  field. 


S 


10 


z  (cm) 

Fig.  2  Fractional  axial  energy  spread,  Yo53z>i»  where  i  =  1,  100,  versus 

distance  along  the  conditioning  waveguide.  The  curves  represent 

3 

100  particles  chosen  randomly  from  a  distribution  of  10  particles. 
In  this  figure,  the  conditioning  field  is  adiabatically  turned  off 


at  =  375  cm. 
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